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Abstract 


The  objective  of  this  project  was  to  validate  the  capability  of  VrWeld  to  simulate  the  weld 
build-up  process  in  two  experimental  setups.  Setup  I  had  a  central  depression  with  dimen¬ 
sions  of  100  x  100  x  3  mm,  while  Setup  II  had  a  central  depression  with  dimensions  of 
200  x  200  x  3  mm. 

This  report  documents  the  validation  of  the  computer  model’s  capability  to  compute  resid¬ 
ual  stress  in  overlay  weld  repairs  of  this  type,  i.e.,  does  the  model  predict  the  measured 
residual  stress  with  sufficient  accuracy,  reliability  and  robustness. 

For  this  project,  the  values  of  the  computed  residual  stress  data  were  mapped  to  the  data 
points  at  which  the  residual  stress  was  measured.  Then  two  types  of  plots  were  made.  The 
first  is  an  image  of  the  stress  field  determined  by  linear  interpolation  onto  the  prefered  grid. 
The  second  is  a  plot  of  curves  of  the  computed  residual  stress  data  and  the  residual  stress 
along  lines  in  the  plane  of  the  panel.  For  both  types,  one  image  of  measured  and  one  image 
of  computed  data  is  shown  on  each  page  to  facilitate  comparison. 

Whether  the  model  predicts  the  measured  residual  stress  with  sufficient  accuracy,  reliability 
and  robustness  for  a  particular  analysis  would  depend  on  the  analysis. 


Resume 


Le  projet  qui  fait  l’objet  du  present  rapport  a  comme  objectif  de  valider  la  capacite  du 
logiciel  VrWeld  de  simuler  le  procede  de  soudage  de  rechargement  de  deux  montages  ex- 
perimentaux.  Le  montage  I  presente  une  depression  centrale  mesurant  100  mm  sur  100  mm 
sur  3  mm,  tandis  que  le  montage  II  en  presente  une  dont  les  dimensions  sont  de  200  mm 
sur  200  mm  sur  3  mm. 

Le  present  rapport  traite  des  travaux  visant  a  valider  la  capacite  du  modele  informatique  de 
calculer  les  contraintes  residuelles  dans  des  structures  reparees  par  ce  type  de  soudage  de 
rechargement ;  en  d’autres  mots,  il  s’agit  d’etablir  si  les  previsions  du  modele,  en  matiere  de 
contraintes  residuelles  mesurees,  sont  adequates  au  chapitre  de  l’exactitude,  de  la  fiabilite 
et  de  la  robustesse. 

Dans  le  cadre  du  present  projet,  on  a  realise  la  mise  en  correspondance  des  valeurs  calculees 
de  contraintes  residuelles  et  des  points  de  donnees  ou  les  contraintes  residuelles  ont  ete 
effectivement  mesurees.  Par  la  suite,  deux  types  de  graphiques  ont  ete  traces.  Ceux  du 
premier  type  constituent  une  image  du  champ  de  contraintes  determine  par  interpolation 
lineaire  sur  le  maillage  privilegie.  Ceux  du  second  type  constituent  un  trace  des  courbes 
des  valeurs  calculees  de  contraintes  residuelles,  le  long  de  lignes  se  trouvant  dans  le  plan 
du  panneau.  Dans  les  deux  cas,  afin  de  faciliter  la  comparaison,  chaque  page  comprend  une 
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image  des  valeurs  mesurees  et  une  autre  des  valeurs  calculees. 

Le  processus  visant  a  determiner  si  le  modele  permet  ou  non  de  prevoir  les  contraintes 
residuelles  mesurees  avec  des  degres  adequats  d’ exactitude,  de  fiabilite  et  de  robustesse, 
est  fonction  de  1’ analyse  realisee. 
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Executive  summary 


Validations  of  Computational  Weld  Models: 

Comparison  of  Residual  Stresses 

John  Goldak;  DRDC  Atlantic  CR  2009-222;  Defence  R&D  Canada  -  Atlantic; 
August  2010. 

Background:  When  welding  is  used  to  to  repair  damage  caused  by  corrosion,  accidents, 
or  re-sealing  interior  accesses  to  the  pressure  hull  of  a  submarine,  changes  in  the  geometry, 
residual  stresses  and  material  properties  have  the  potential  to  affect  its  in-service  behaviour. 
In  particular,  there  is  a  need  to  develop  the  capability  to  estimate  the  extent  to  which  weld¬ 
ing  changes  the  resistance  of  the  pressure  hull  to  buckling,  fatigue  and  fracture. 

To  this  end,  DRDC  Atlantic  commissioned  BMT  Fleet  Technologies  Limited  to  make  over¬ 
lay  welds  on  two  rigid  panels.  The  structures  were  instrumented  with  thermocouples  and 
strain  gauges  which  recorded  the  change  in  temperature  and  strains  during  the  preheating, 
welding,  and  cool-down.  The  distortion  of  the  bottom  of  the  plate  in  the  final  structure  was 
measured  and  reported  in  Monitoring  of  Welding  Induced  Strains  during  Cladding  [1]. 

DRDC  Atlantic  commissioned  Goldak  Technologies  Inc.  to  model  the  overlay  welding 
procedure  used  to  make  the  two  panels  to  predict  the  transient  distortion,  strain  and  stress 
in  the  two  rigid  panels.  These  results  are  summarized  in  the  report  Thermo-mechanical 
modelling  of  welding  induced  strains:  Numerical  validation  of  the  weld  build-up  process 
[2]. 

This  document  reports  on  the  comparison  of  residual  stresses  in  the  two  panels  measured  by 
neutron  diffraction  and  computed  by  the  computational  weld  mechanics  software  VrWelcl. 

Principal  results: 

1.  The  100  x  100  mm  built-up  region  were  simulated  using  27  overlay  weld  passes 
while  the  200  x  200  mm  built-up  region  was  simulated  with  51  overlay  weld  passes. 

2.  The  influence  of  filler  metal  chemistry  was  also  included  for  two  cases. 

3.  The  residual  stresses  in  experimental  setup  I  with  the  100  x  100  mm  built-up  region 
and  experimental  setup  II  with  the  200  x  200  mm  built-up  region  computed  by  a 
computational  weld  mechanics  model,  VrWeld,  and  measured  by  neutron  diffraction 
were  compared.  In  the  experimetal  setup  the  resiidual  stresses  were  were  measured 
at  points  in  the  longitudinal  and  transverse  planes. 
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4.  While  a  visual  examination  of  the  residual  stress  plots  suggest  that  there  is  some 
agreement  between  the  measured  and  computed  values  of  residual  stress,  to  the  au¬ 
thor’s  knowledge,  no  generally  accepted  methods  to  compute  a  quantitative  measure 
of  the  agreement  are  available. 

Significance  of  results:  To  validate  the  computer  model,  an  application  must  be  defined 
and  the  sensitivity  of  the  model  results  to  that  application  must  be  determined.  Whether  this 
agreement  is  sufficient  for  a  particular  design  decision  would  depend  on  the  applications. 

Future  work:  Based  on  past  and  current  experience,  the  following  recommendation  and 
future  work  are  suggested: 

A  virtual  Design  of  Experiment  (DOE)  should  be  done  to  determine  the  sensitivity  of  the 
solution  to  model  parameters.  This  would  enable  the  agreement  between  measured  and 
computed  residual  stress  to  be  improved.  The  computed  longitudinal  stresses  in  the  weld 
metal  appear  to  be  higher  than  the  measured  longitudinal  stresses.  This  higher  tensile 
residual  stress  in  the  weld  metal  could  generate  the  higher  longitudinal  compressive  stress 
just  below  the  high  tensile  stress.  This  suggests  that  it  would  be  worth  doing  a  DOE  analysis 
to  compute  the  sensitivity  to  changes  in  the  yield  stress  of  martensite  and  bainite  phases  to 
determine  if  this  improves  the  agreement  between  measured  and  computed  residual  stress. 

The  effect  of  martensite  on  the  tempering  to  reduce  the  hardness  of  martensite  formed  in 
one  weld  pass  by  a  later  pass  should  be  analyzed.  This  reduction  in  hardness  suggests  that 
the  yield  stress  of  martensite  be  reduced  from  approximately  900  MPa  of  the  untempered 
martensite.  If  a  model  is  developed  for  tempering  of  martensite,  then  the  model  could  be 
included  in  the  computer  model. 

Hardness  measurements  of  the  plate  suggest  that  the  hardness  has  a  significant  fluctuations 
in  space,  i.e.,  it  is  not  homogeneous.  It  would  be  of  interest  to  study  the  effect  of  such 
fluctuations  on  the  residual  stress.  Transformation  plasticity  is  a  result  of  fluctuations  at  the 
micro-scale.  However,  the  author  is  not  aware  of  investigations  of  the  effects  of  fluctuations 
at  the  macro-scale  level. 

The  meshing  could  be  further  improved,  particularly  the  stopping  points  of  each  weld  pass. 
Goldak  Technologies  Inc.  is  exploring  ways  to  do  this,  such  as  domain  decomposition  with 
a  small  domain  with  a  fine  mesh  and  space-time  FEM  analysis  with  a  graded  mesh. 

Difficulties  Encountered:  It  would  be  useful  to  establish  a  reference  coordinate  system 
and  then  refer  all  experimental  and  computed  data  to  the  reference  coordinate  system. 
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Validations  of  Computational  Weld  Models: 

Comparison  of  Residual  Stresses 

John  Goldak;  DRDC  Atlantic  CR  2009-222;  R  &  D  pour  la  defense  Canada  - 
Atlantique;  aout2010. 

Introduction:  Lorsque  des  travaux  de  soudage  sont  realises  afin  de  reparer  des  dommages 
causes  par  la  corrosion  ou  des  accidents,  ou  pour  etancheiser  de  nouveau  les  acces  in- 
terieurs  de  la  coque  epaisse  d’un  sous-marin,  le  comportement  en  service  des  structures 
peut  subir  les  effets  de  la  modification  de  la  geometrie,  de  la  nature  et  de  rimportance  des 
contraintes  residuelles  et  des  proprietes  des  materiaux.  II  existe  done  un  besoin  particulier 
dans  ce  domaine,  soit  l’elaboration  d’un  outil  ayant  la  capacite  d’estimer  l’ampleur  des 
changements  que  provoque  le  soudage  au  chapitre  de  la  resistance  d’une  coque  epaisse  au 
gauchissement,  a  la  fatigue  et  a  la  rupture. 

C’est  pourquoi  RDDC  Atlantique  a  retenu  les  services  de  la  societe  BMT  Fleet  Techno¬ 
logies  Limited  afin  que  celle-ci  execute  le  soudage  de  rechargement  de  deux  panneaux 
rigides.  Les  structures  etaient  munies  de  thermocouples  et  de  jauges  de  deformation  per- 
mettant  d’enregistrer  les  variations  de  temperature  et  les  deformations  durant  les  etapes  de 
prechauffage,  de  soudage  et  de  refroidissement.  La  distorsion  de  la  partie  inferieure  de  la 
plaque  soudee  finale  a  ete  mesuree  et  les  resultats  sont  presentes  dans  le  rapport  intitule 
Monitoring  of  Welding  Induced  Strains  during  Cladding  [1], 

RDDC  Atlantique  a  aussi  retenu  les  services  de  la  societe  Goldak  Technologies  Inc.  afin 
que  celle-ci  execute  la  modelisation  du  procede  de  soudage  de  rechargement  utilise  pour 
fabriquer  les  deux  panneaux  rigides.  Un  resume  des  resultats  apparait  dans  le  rapport 
Thermo-mechanical  modelling  of  welding  induced  strains  :  Numerical  validation  of  the 
weld  build-up  process  [2]. 

Le  present  document  traite  de  la  comparaison  des  contraintes  residuelles  que  presentent  les 
deux  panneaux,  mesurees  par  diffraction  des  neutrons,  et  de  celles  calculees  au  moyen  du 
logiciel  VrWeld,  qui  simule  la  mecanique  des  procedes  de  soudage. 

Resultats: 

1.  La  simulation  de  la  zone  de  rechargement  mesurant  100  mm  sur  100  mm  a  ete  reali- 
see  en  employant  27  passes  de  soudure  de  rechargement,  tandis  que  celle  de  la  zone 
mesurant  200  mm  sur  200  mm  a  ete  executee  en  employant  51  passes. 

2.  On  a  aussi  tenu  compte  de  la  composition  chimique  du  metal  d’apport,  et  ce,  dans 
les  deux  cas  etudies. 
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3.  On  a  compare  les  valeurs  des  contraintes  residuelles  du  montage  experimental  I, 
presentant  une  zone  de  rechargement  de  100  mm  sur  100  mm,  et  du  montage  expe¬ 
rimental  II,  presentant  une  zone  de  rechargement  de  200  mm  sur  200  mm,  calculees 
au  moyen  du  logiciel  de  modelisation  de  la  mecanique  des  precedes  de  soudage 
VrWeld,  et  celles  mesurees  par  diffraction  des  neutrons.  Les  contraintes  residuelles 
mesurees  dans  les  montages  experimentaux  l’ont  ete  en  des  points  particuliers  des 
plans  longitudinal  et  transversal. 

4.  Un  examen  visuel  des  graphiques  des  contraintes  residuelles  laisse  supposer  qu’il 
existe  un  certain  accord  entre  les  valeurs  de  contraintes  residuelles  mesurees  et  cal¬ 
culees,  mais  a  la  connaissance  de  1’ auteur  du  present  rapport,  il  n’ existe  cependant 
aucune  methode  eprouvee  qui  permet  d’evaluer  quantitativement  le  degre  de  cet  ac¬ 
cord  ou  de  la  conformite  des  differents  resultats. 

Portee:  Afin  d’effectuer  la  validation  du  modele  numerique,  il  est  essentiel  de  definir  une 
application  et  de  determiner  la  sensibilite  des  resultats  de  modelisation  par  rapport  a  cette 
application.  Le  processus  visant  a  determiner  si  l’accord  obtenu  est  suffisant  pour  prendre 
des  decisions  de  conception  particulieres  serait  fonction  des  applications  envisagees. 

Recherches  futures:  D’apres  l’experience  acquise  dans  le  cadre  du  present  projet  et  de 
projets  anterieurs,  on  peut  recommander  l’execution  des  travaux  et  activites  suivantes  : 

Un  dispositif  experimental  (DE)  virtuel  devrait  etre  elabore  afin  de  determiner  la  sensibilite 
de  la  solution  par  rapport  aux  parametres  du  modele.  Cette  approche  permettrait  d’amelio- 
rer  l’accord  entre  les  valeurs  de  contraintes  residuelles  mesurees  et  calculees.  Les  valeurs 
calculees  des  contraintes  longitudinales  dans  le  metal  d’apport  (metal  de  soudure)  semblent 
etre  superieures  a  celles  mesurees.  Ces  contraintes  residuelles  de  traction  plus  elevees  dans 
le  metal  d’apport  pourraient  constituer  la  cause  des  contraintes  residuelles  de  compression 
plus  elevees  observees  dans  la  zone  inferieure  adjacente.  Ces  faits  laissent  supposer  qu’il 
serait  utile  de  realiser  une  analyse  de  DE  afin  de  calculer  la  sensibilite  aux  changements 
de  la  contrainte  d’ecoulement  dans  les  phases  de  martensite  et  de  bainite  et  de  determi¬ 
ner  s’il  est  ainsi  possible  d’ameliorer  l’accord  entre  les  valeurs  mesurees  et  calculees  des 
contraintes  residuelles. 

Il  faudrait  analyser  les  effets  de  la  martensite  sur  le  processus  de  recuit  afin  de  reduire  la 
durete  de  la  phase  martensitique  qui  est  formee  dans  une  passe  de  soudure  par  une  passe 
ulterieure.  Une  telle  reduction  de  la  durete  semble  indiquer  que  la  contrainte  d’ecoulement 
de  la  martensite  pourrait  etre  reduite,  par  rapport  a  une  valeur  de  quelque  900  MPa  pour  la 
martensite  n’ayant  pas  subi  un  recuit.  Si  un  modele  du  processus  de  recuit  de  la  martensite 
peut  etre  elabore,  il  pourrait  etre  integre  au  modele  informatique. 

Les  mesures  de  la  durete  de  la  plaque  laissent  supposer  que  cette  propriete  varie  grandement 
dans  l’espace,  c.-a-d.  qu’elle  est  inhomogene.  Il  serait  interessant  d’etudier  les  effets  de 
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ces  variations  sur  les  contraintes  residuelles.  La  plastic  ite  de  transformation  decoule  de 
variations  qui  se  produisent  a  l’echelle  du  micrometre.  A  la  connaissance  de  l’auteur  du 
present  rapport,  il  n’existe  cependant  aucune  etude  portant  sur  les  effets  des  variations  a 
l’echelle  macro scopique. 

La  qualite  du  maillage  pourrait  etre  amelioree,  particulierement  aux  points  d’ interruption 
de  chaque  passe  de  soudure.  La  societe  Goldak  Technologies  Inc.  examine  actuellement 
des  moyens  d’atteindre  cet  objectif,  par  exemple  en  employant  un  domaine  de  petites  di¬ 
mensions  a  maillage  fin  et  T  analyse  FEM  (modele  par  elements  finis)  dans  l’espace  et  le 
temps  avec  maillage  dimensionne. 

Problemes  rencontres:  II  serait  utile  d’etablir  un  systeme  de  coordonnees  de  reference  et 
d’exprimer  toutes  les  donnees  experimentales  et  calculees  en  fonction  de  ce  systeme. 
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Figure  1 1 :  The  images  are  for  longitudinal  stress  on  the  small  panel  in  the  plane 
parallel  to  the  welding  direction.  In  the  top  image  the  FEM  data  is 
plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 
the  FEM  data .  11 

Figure  12:  The  images  are  for  transverse  stress  on  the  small  panel  in  the  plane 
parallel  to  the  welding  direction.  In  the  top  image  the  FEM  data  is 
plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 
the  FEM  data .  12 

Figure  13:  The  images  are  for  normal  stress  on  the  small  panel  in  the  plane  parallel 
to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with 
linear  interpolation  between  data  points  on  a  cut-plane  mesh  in  the  3D 
mesh  used  to  solve  for  stresses.  In  the  remaining  images  data  is  plotted 
using  linear  interpolation  between  data  points  on  the  grid  of  points  used 
to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks. 

The  second  image  from  top  shows  the  measured  data.  The  third  image 

from  the  top  shows  the  FEM  data  mapped  to  the  grid.  The  bottom 

image  shows  the  difference  of  the  measured  data  from  the  FEM  data.  .  .  13 

Figure  14:  The  images  are  for  longitudinal  stress  on  the  small  panel  in  the  plane 
perpendicular  to  the  welding  direction.  In  the  top  image  the  FEM  data 
is  plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 


the  FEM  data .  14 
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Figure  15:  The  images  are  for  transverse  stress  on  the  small  panel  in  the  plane 

perpendicular  to  the  welding  direction.  In  the  top  image  the  FEM  data 
is  plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 
the  FEM  data .  15 

Figure  16:  The  images  are  for  normal  stress  on  the  small  panel  in  the  plane 

perpendicular  to  the  welding  direction.  In  the  top  image  the  FEM  data 
is  plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 
the  FEM  data .  16 

Figure  17:  The  images  are  for  longitudinal  stress  on  the  large  panel  in  the  plane 
parallel  to  the  welding  direction.  In  the  top  image  the  FEM  data  is 
plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 
the  FEM  data .  17 

Figure  18:  The  images  are  for  transverse  stress  on  the  large  panel  in  the  plane 
parallel  to  the  welding  direction.  In  the  top  image  the  FEM  data  is 
plotted  with  linear  interpolation  between  data  points  on  a  cut-plane 
mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the  remaining 
images  data  is  plotted  using  linear  interpolation  between  data  points  on 
the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are 
shown  as  black  disks.  The  second  image  from  top  shows  the  measured 
data.  The  third  image  from  the  top  shows  the  FEM  data  mapped  to  the 
grid.  The  bottom  image  shows  the  difference  of  the  measured  data  from 
the  FEM  data .  18 
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Figure  19:  The  images  are  for  normal  stress  on  the  large  panel  in  the  plane  parallel 
to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with 
linear  interpolation  between  data  points  on  a  cut-plane  mesh  in  the  3D 
mesh  used  to  solve  for  stresses.  In  the  remaining  images  data  is  plotted 
using  linear  interpolation  between  data  points  on  the  grid  of  points  used 
to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks. 

The  second  image  from  top  shows  the  measured  data.  The  third  image 

from  the  top  shows  the  FEM  data  mapped  to  the  grid.  The  bottom 

image  shows  the  difference  of  the  measured  data  from  the  FEM  data.  .  .  19 

Figure  20:  Measured  (top)  and  computed  (bottom)  longitudinal  residual  stress  in 

the  small  panel  in  the  plane  parallel  to  the  weld  direction  is  shown.  ...  20 

Figure  21:  Measured  (top)  and  computed  (bottom)  transverse  residual  stress  in  the 

small  panel  in  the  plane  parallel  to  the  weld  direction  is  shown . 21 

Figure  22:  Measured  (top)  and  computed  (bottom)  normal  residual  stress  in  the 

small  panel  in  the  plane  parallel  to  the  weld  direction  is  shown . 22 

Figure  23:  Measured  (top)  and  computed  (bottom)  longitudinal  residual  stress  in 

the  small  panel  in  the  plane  normal  to  the  weld  direction  is  shown.  ...  23 

Figure  24:  Measured  (top)  and  computed  (bottom)  transverse  residual  stress  in  the 

small  panel  in  the  plane  normal  to  the  weld  direction  is  shown . 24 

Figure  25:  Measured  (top)  and  computed  (bottom)  normal  residual  stress  in  the 

small  panel  in  the  plane  normal  to  the  weld  direction  is  shown . 25 

Figure  26:  Measured  (top)  and  computed  (bottom)  longitudinal  residual  stress  in 

the  large  panel  in  the  plane  parallel  to  the  weld  direction  is  shown.  ...  26 

Figure  27:  Measured  (top)  and  computed  (bottom)  transverse  residual  stress  in  the 

large  panel  in  the  plane  parallel  to  the  weld  direction  is  shown . 27 

Figure  28:  Measured  (top)  and  computed  (bottom)  normal  residual  stress  in  the 

large  panel  in  the  plane  parallel  to  the  weld  direction  is  shown . 28 

Figure  29:  The  difference  between  the  computed  and  measured  transverse  residual 
stress  in  the  large  panel  in  the  plane  parallel  to  the  weld  direction  is 
shown.  Compare  to  Fig.  27 . 30 
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1  Modelling  the  Structures 


The  fine  mesh  of  the  100  x  100  overlay  has  15,606  elements  and  24,727  nodes.  The  fine 
mesh  of  the  200  x  200  overlay  has  50,808  elements  and  69,627  nodes.  The  coarse  mesh 
of  the  200  x  200  overlay  has  21,978  elements  and  32,880  nodes.  The  fine  meshes  are 
considered  to  be  very  high  quality  meshes.  The  cross-section  of  each  filler  metal  weld  pass 
had  from  6  to  15  brick  elements.  The  mesh  for  the  region  being  welded  including  the  filler 
metal  were  meshed  with  8-node  bricks.  The  complement  of  the  region  being  welded  was 
meshed  with  6-node  prism  and  8-node  brick  elements.  The  meshing  was  done  in  VrSuite. 
Various  views  of  these  meshes  are  shown  in  Figures  1-5. 


Figure  1:  A  cross-section  of  the  fine  mesh  of  the  100  x  100  overlay  is  shown. 


Figure  2:  A  cross-section  of  the  fine  mesh  of  the  100  x  100  overlay  is  shown.  The  filler 
metal  cross-sections  are  coloured. 
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Figure  3:  A  cross-section  of  the  fine  mesh  of  the  200  x  200  overlay  is  shown. 


Figure  4:  A  cross-section  of  the  fine  mesh  of  the  200  x  200  overlay  is  shown. 
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Figure  5:  A  view  of  the  fine  mesh  of  the  100  x  100  overlay  structure  is  shown.  The  green 
ball  represents  a  virtual  thermocouple.  The  pairs  of  white  balls  represent  virtual  strain 
gauges. 
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1.1  Thermal  Analysis 

Weld  passes  were  deposited  in  the  depressions  in  a  symmetrical  sequence  which  incremen¬ 
tally  increased  the  clad  width  and  hence  balances  the  welding  induced  strain.  The  plate 
material  was  HY80  and  the  filler  metal  Fortrex  AWS  E9016.  During  the  welding  process 
transient  temperatures  were  recorded  with  thermocouples  and  transient  strains  with  strain 
gauges.  After  welding  and  cool-down  to  room  temperature,  the  distortion  on  the  bottom  of 
each  plate  was  measured  along  with  through-thickness  residual  stress  measurements.  CR 
2008-283  [2]  provides  further  details  of  the  analysis  procedure. 

1.2  Chemical  Compositions 

Two  variations  in  the  chemical  compositions  were  examined.  In  the  first  case  both  the  base 
metal  and  weld  filler  metal  utilized  the  parent  metal  composition  listed  in  Table  1,  while  in 
the  second  case,  the  weld  metal  utilized  the  as-deposited  chemical  composition.  . 

Table  1:  Chemical  composition  of  parent  and  consumable  materials  (in  wt%) 


Composition 

C 

Mn 

Si 

S 

P 

Ni 

Cr 

Mo 

Parent  metal 

0.14 

0.27 

0.15 

0.002 

0.006 

2.2 

1.1 

0.25 

Weld  metal  as-deposited 

0.08 

1.3 

0.17 

0.002 

0.009 

1.1 

0.23 

0.18 
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2  Experimental  Measurements 


McLaughlin  [8]  obtained  neutron  diffraction  measurements  from  the  weld  built-up  panels 
along  the  four  planes  illustrated  in  Figure  6  and  7  for  the  small  and  large  weld  patches, 
respectively.  The  experimental  residual  stresses  were  measured  on  a  6  x  24  grid  on  a  lon¬ 
gitudinal  plane  and  a  7  x  13  grid  on  transverse  plane  on  the  small  weld  patch  and  a  7  x  13 
grid  on  a  longitudinal  plane  on  the  large  weld  patch.  These  measurement  locations  are 
represented  by  red  dots  in  Figures  6  and  7. 


Figure  6:  The  small  experimental  calibration  panel  showing  the  approximate  location  of 
the  two  through  thickness  planes  used  to  acquire  the  residual  stress  measurements.  Along 
the  longitudinally  oriented  plane,  the  measurements  extended  across  the  whole  weld  panel, 
while  for  the  transverse  plane,  only  half  of  the  build-up  panel  was  measured.  Measurements 
were  acquired  in  the  longitudinal,  transverse  and  normal  direction  at  geometrically  similar 
positions. 
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Figure  7:  The  large  experimental  calibration  panel  showing  the  approximate  location  of 
the  two  through  thickness  planes  used  to  acquire  the  residual  stress  measurements.  Along 
the  longitudinally  oriented  plane,  the  measurements  extended  across  the  whole  weld  panel, 
while  for  the  transverse  plane,  only  half  of  the  build-up  panel  was  measured. 


6 


DRDC  Atlantic  CR  2009-222 


3  Influence  of  Weld  Composition 


The  influence  that  the  weld  metal  composition  has  on  the  computed  residual  stresses  are 
examined  in  Figure  8  to  10.  These  compare  the  longitudinal,  transverse  and  normal  stress 
components  along  a  transverse  weld  section  on  the  100  x  100  weld  patch  (See  Figure  6). 
In  all  cases  subtle  changes  in  the  residual  stress  field  are  calculate  demonstrating  that  the 
influence  of  chemical  composition  on  mechanical  properties  can  be  computed. 
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(a)  Weld  uses  base  composition 


StressNodal(MPa) 


(b)  As-deposited  composition 

Figure  8:  Influence  of  the  weld  metal  chemistry  on  the  computed  residual  stresses  (100  x 
100  Longitudinal  Stress  Component  Across  the  Weld) 
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StrenNodatCMPa) 


(a)  Weld  uses  base  composition 


StressNodal(MPa) 


(b)  As-deposited  composition 

Figure  9:  Influence  of  the  weld  metal  chemistry  on  the  computed  residual  stresses  (100  x 
100  Transverse  Stress  Component  Across  the  Weld) 
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(a)  Weld  uses  base  composition 
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(b)  As-deposited  composition 

Figure  10:  Influence  of  the  weld  metal  chemistry  on  the  computed  residual  stresses  (100  x 
100  Normal  Stress  Component  Across  the  Weld) 
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4  Comparison  of  Measured  and  Computed 
Residual  stress 


The  longitudinal,  transverse  and  normal  residual  stress  is  measured  at  points  on  a  coarse 
grid  in  a  longitudinal  and  a  transverse  plane  in  the  small  patch  (Figure  6)  and  in  a  longitu¬ 
dinal  plane  in  the  large  patch  (Figure  7).  The  images  of  the  measured  stress  distribution  are 
shown  in  the  second  image  from  the  top  in  Figures  1 1  to  19.  These  images  were  computed 
by  interpolating  data  from  this  coarse  grid  to  the  the  pixels  in  the  image. 

The  longitudinal,  transverse  and  normal  residual  stress  is  computed  at  points  in  a  much 
finer  FEM  mesh.  The  images  of  the  stress  distribution  are  computed  by  interpolating  data 
from  this  fine  FEM  mesh  to  the  the  pixels  in  the  top  image  on  Figures  1 1  to  19.  Note  these 
images  do  not  use  interpolation  on  the  coarse  grid  to  the  pixels  in  the  image. 

An  FEM  mesh  was  created  in  which  the  nodes  correspond  to  data  points  used  to  measure 
the  residual  stress.  Data  from  the  fine  FEM  mesh  was  mapped  to  these  coarse  grid  points 
and  then  interpolated  to  image  pixels  to  create  the  third  from  the  top  image  in  on  Figures 
11  to  19  .  These  images  were  computed  by  interpolating  data  from  this  coarse  grid  to  the 
the  pixels  in  the  image. 

The  bottom  image  in  Figures  11  to  19  are  the  difference  between  the  measured  and  com¬ 
puted  stress  at  the  grid  points  which  is  interpolated  to  pixels  to  create  the  image.  These 
images  were  computed  by  interpolating  data  from  this  coarse  grid  to  the  the  pixels  in  the 
image. 

Figures  20  to  25  show  line  plots  for  measured  and  computed  stresses  on  each  line  of  the 
associated  coarse  grid  in  the  small  panel.  Figures  26  to  28  show  line  plots  for  measured 
and  computed  stresses  on  each  line  of  the  associated  coarse  grid  in  the  large  panel. 
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Figure  11:  The  images  are  for  longitudinal  stress  on  the  small  panel  in  the  plane  parallel 
to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  12:  The  images  are  for  transverse  stress  on  the  small  panel  in  the  plane  parallel  to 
the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  13:  The  images  are  for  normal  stress  on  the  small  panel  in  the  plane  parallel  to 
the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  14:  The  images  are  for  longitudinal  stress  on  the  small  panel  in  the  plane  per¬ 
pendicular  to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear 
interpolation  between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for 
stresses.  In  the  remaining  images  data  is  plotted  using  linear  interpolation  between  data 
points  on  the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as 
black  disks.  The  second  image  from  top  shows  the  measured  data.  The  third  image  from 
the  top  shows  the  FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference 
of  the  measured  data  from  the  FEM  data. 
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Figure  15:  The  images  are  for  transverse  stress  on  the  small  panel  in  the  plane  perpen¬ 
dicular  to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear 
interpolation  between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for 
stresses.  In  the  remaining  images  data  is  plotted  using  linear  interpolation  between  data 
points  on  the  grid  of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as 
black  disks.  The  second  image  from  top  shows  the  measured  data.  The  third  image  from 
the  top  shows  the  FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference 
of  the  measured  data  from  the  FEM  data. 
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Figure  16:  The  images  are  for  normal  stress  on  the  small  panel  in  the  plane  perpendicular 
to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  17:  The  images  are  for  longitudinal  stress  on  the  large  panel  in  the  plane  parallel 
to  the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  18:  The  images  are  for  transverse  stress  on  the  large  panel  in  the  plane  parallel  to 
the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  19:  The  images  are  for  normal  stress  on  the  large  panel  in  the  plane  parallel  to 
the  welding  direction.  In  the  top  image  the  FEM  data  is  plotted  with  linear  interpolation 
between  data  points  on  a  cut-plane  mesh  in  the  3D  mesh  used  to  solve  for  stresses.  In  the 
remaining  images  data  is  plotted  using  linear  interpolation  between  data  points  on  the  grid 
of  points  used  to  measure  residual  stresses.  The  grid  points  are  shown  as  black  disks.  The 
second  image  from  top  shows  the  measured  data.  The  third  image  from  the  top  shows  the 
FEM  data  mapped  to  the  grid.  The  bottom  image  shows  the  difference  of  the  measured 
data  from  the  FEM  data. 
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Figure  20:  Measured  (top)  and  computed  (bottom)  longitudinal  residual  stress  in  the  small 
panel  in  the  plane  parallel  to  the  weld  direction  is  shown. 
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Figure  21:  Measured  (top)  and  computed  (bottom)  transverse  residual  stress  in  the  small 
panel  in  the  plane  parallel  to  the  weld  direction  is  shown. 
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Figure  22:  Measured  (top)  and  computed  (bottom)  normal  residual  stress  in  the  small 
panel  in  the  plane  parallel  to  the  weld  direction  is  shown. 
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Figure  23:  Measured  (top)  and  computed  (bottom)  longitudinal  residual  stress  in  the  small 
panel  in  the  plane  normal  to  the  weld  direction  is  shown. 
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Figure  24:  Measured  (top)  and  computed  (bottom)  transverse  residual  stress  in  the  small 
panel  in  the  plane  normal  to  the  weld  direction  is  shown. 
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Figure  25:  Measured  (top)  and  computed  (bottom)  normal  residual  stress  in  the  small 
panel  in  the  plane  normal  to  the  weld  direction  is  shown. 
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Figure  26:  Measured  (top)  and  computed  (bottom)  longitudinal  residual  stress  in  the  large 
panel  in  the  plane  parallel  to  the  weld  direction  is  shown. 
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Figure  27:  Measured  (top)  and  computed  (bottom)  transverse  residual  stress  in  the  large 
panel  in  the  plane  parallel  to  the  weld  direction  is  shown. 
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Figure  28:  Measured  (top )  and  computed  (bottom )  normal  residual  stress  in  the  large  panel 
in  the  plane  parallel  to  the  weld  direction  is  shown. 
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5  Discussion  of  Comparison  of  Measured 
and  Computed  Residual  Stress 


The  validation  of  the  code  for  a  computer  model  requires  comparing  computed  results  with 
experimental  results.  However,  no  rules  have  been  established  that  define  the  best  methods 
for  measuring  the  difference  between  measured  and  computed  results.  The  situation  is 
particularly  challenging  for  3D  continuum  fields  such  as  the  field  for  the  residual  stress 
tensor.  Sensors  sample  a  finite  number  of  points  and  each  point  integrates  over  a  finite 
volume  of  space.  The  computer  model  must  discretize  space  and  time.  From  the  sensor  data 
points  and  the  computer  model  data,  one  would  like  to  compute  3D  residual  stress  fields 
and  evaluate  the  ‘distance’  between  the  measured  and  computed  data  fields.  In  mathematics 
the  notion  of  distance  between  two  functions  or  fields  requires  a  norm.  There  is  little 
or  no  agreement  in  the  published  literature  on  which  norm  or  norms  are  best.  It  is  the 
responsibility  of  the  decision  maker  to  decide  if  the  agreement  between  the  measured  data 
and  data  computed  by  the  model  is  sufficient  to  justify  use  of  the  model  for  a  particular 
application  [6]  [9] . 

The  numerical  methodology  for  comparing  measured  and  computed  values  of  a  3D  con¬ 
tinuum  field  such  as  a  stress  tensor  is  a  relatively  immature.  Figure  27  shows  curves  for 
computed  and  measured  transverse  stress.  A  person  looking  at  the  curves,  might  say  the 
agreement  is  quite  good.  Figure  28  shows  the  difference  between  measured  and  computed 
transverse  stresses  that  are  plotted  in  Figure  27.  The  maximum  differences  range  from 
-250MPa  to  270  MPa  while  the  stresses  in  Figure  27  range  from  -350  to  +400  MPa. 

For  this  project,  the  values  of  the  computed  residual  stress  data  were  mapped  to  the  data 
points  at  which  the  residual  stress  was  measured.  Then  two  types  of  images  were  made. 
One  image  type  is  a  jpg  of  the  stress  field  determined  by  linear  interpolation  on  the  grid 
defined  by  the  measured  data  points.  The  second  image  type  is  a  plot  of  curves  of  the 
computed  residual  stress  data  and  the  residual  stress  along  lines  in  the  plane  of  the  panel. 
For  both  types,  one  image  of  measured  and  one  image  of  computed  data  is  shown  on  each 
page  to  facilitate  comparison.  The  residual  stress  measured  at  points  on  a  coarse  grid  has 
a  number  of  uncertainties  associated  with  the  neutron  diffraction  technique,  the  welding 
procedure  and  the  manufacture  of  the  plate.  The  computed  values  of  the  residual  stress 
also  have  a  number  of  uncertainties  associated  with  the  analysis.  In  the  author’s  opinion, 
the  greatest  uncertainties  in  the  computed  values  of  the  residual  stress  are  associated  with 
the  material  properties  that  are  functions  of  temperature,  microstructure  and  time.  There  is 
also  a  significant  discretization  error  caused  by  using  a  finite  mesh. 

If  the  measured  longitudinal,  transverse  and  normal  residual  stresses  are  assumed  to  be 
principal  stresses  and  the  associated  Cartesian  directions  are  assumed  to  be  the  associated 
eigenvectors,  then  the  measured  values  define  the  field  of  the  stress  tensor.  The  com¬ 
puted  residual  stress  does  not  assume  that  the  longitudinal,  transverse  and  normal  residual 
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Figure  29:  The  difference  between  the  computed  and  measured  transverse  residual  stress 
in  the  large  panel  in  the  plane  parallel  to  the  weld  direction  is  shown.  Compare  to  Fig.  21. 


stresses  are  principal  stresses.  See  [2]  for  images  of  the  computed  principal  stresses.  It 
would  be  possible  to  compute  the  deviation  of  the  computed  principal  stresses  from  the 
Cartesian  axes  and  thus  assess  the  error  in  assuming  the  Cartesian  directions  are  directions 
of  principal  stresses. 

The  computed  stresses  are  kinematically  admissible  stresses,  i.e.,  the  displacement  field  is 
C°  continuous,  satisfies  Dirichlet  or  prescribed  boundary  conditions,  the  strain  is  a  func¬ 
tion  of  the  symmetrized  displacement  gradient  and  the  4th  order  elasticity-plasticity  tensor 
maps  the  strain  tensor  to  the  stress  tensor.  A  statically  admissible  stress  field  would  be  C° 
continuous,  satisfies  Neumann  or  traction  boundary  conditions  and  the  momentum  equa¬ 
tion.  In  particular,  a  statically  admissible  stress  field  is  in  equilibrium  with  zero  traction 
boundary  conditions.  The  Prager-Synge  hyper-circle  theorem  states  that  the  exact  solution 
is  bounded  between  the  kinematically  admissible  and  the  statically  admissible  solutions. 
See  Mashie  [7]  for  further  details. 
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6  Recommendations  and  Future  Work 


Based  on  past  and  current  experience,  the  following  recommendation  and  future  work  are 
suggested: 

A  virtual  Design  of  Experiment  (DOE)  should  be  done  to  determine  the  sensitivity  of  the 
solution  to  model  parameters.  This  would  enable  the  agreement  between  measured  and 
computed  residual  stress  to  be  improved.  The  computed  longitudinal  stresses  in  the  weld 
metal  appear  to  be  higher  than  the  measured  longitudinal  stresses.  This  higher  tensile 
residual  stress  in  the  weld  metal  could  generate  the  higher  longitudinal  compressive  stress 
just  below  the  high  tensile  stress.  This  suggests  that  it  would  be  worth  doing  a  DOE  analysis 
to  compute  the  sensitivity  to  changes  in  the  yield  stress  of  martensite  and  bainite  phases  to 
determine  if  this  improves  the  agreement  between  measured  and  computed  residual  stress. 

The  effect  of  martensite  on  the  tempering  to  reduce  the  hardness  of  martensite  formed  in 
one  weld  pass  by  a  later  pass  should  be  analyzed.  This  reduction  in  hardness  suggests  that 
the  yield  stress  of  martensite  be  reduced  from  approximately  900  MPa  of  the  untempered 
martensite.  If  a  model  is  developed  for  tempering  of  martensite,  then  the  model  could  be 
included  in  the  computer  model. 

Hardness  measurements  of  the  plate  suggest  that  the  hardness  has  a  significant  fluctuations 
in  space,  i.e.,  it  is  not  homogeneous.  It  would  be  of  interest  to  study  the  effect  of  such 
fluctuations  on  the  residual  stress.  Transformation  plasticity  is  a  result  of  fluctuations  at  the 
micro-scale.  However,  the  author  is  not  aware  of  investigations  of  the  effects  of  fluctuations 
at  the  macro-scale  level. 

The  meshing  could  be  further  improved,  particularly  the  stopping  points  of  each  weld  pass. 
Goldak  Technologies  Inc.  is  exploring  ways  to  do  this,  such  as  domain  decomposition  with 
a  small  domain  with  a  fine  mesh  and  space-time  FEM  analysis  with  a  graded  mesh. 
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